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Abstract 

Background: Protein pliospliatases (PPs) play critical roles in various cellular processes through the reversible 
protein phosphorylation that dictates many signal transduction pathways among organisms. Recently, PPs in 
Arabidopsis and rice have been identified, while the whole complement of PPs in maize is yet to be reported. 

Results: In this study, we have identified 159 PP-encoding genes in the maize genome. Phylogenetic analyses 
categorized the ZmPP gene family into 3 classes (PP2C, PTP, and PP2A) with considerable conservation among 
classes. Similar intron/exon structural patterns were observed in the same classes. Moreover, detailed gene 
structures and duplicative events were then researched. The expression profiles of Zm PPs under different 
developmental stages and abiotic stresses (including salt, drought, and cold) were analyzed using microarray and 
RNA-seq data. A total of 152 members were detected in 18 different tissues representing distinct stages of maize 
plant developments. Under salt stress, one gene was significantly up-expressed in seed root (SR) and one gene was 
down-expressed in primary root (PR) and crown root (CR), respectively. As for drought stress condition, 13 genes 
were found to be differentially expressed in leaf out of which 10 were up-regulated and 3 exhibited down-regulation. 
Additionally, 13 up-regulated and 3 down-regulated genes were found in cold-tolerant line ETH-DH7. Furthermore, 
real-time PCR was used to confirm the expression patterns of ZmPPs. 

Conclusions: Our results provide new insights into the phylogenetic relationships and characteristic functions of 
maize PPs and will be useful in studies aimed at revealing the global regulatory network in maize abiotic stress 
responses, thereby contributing to the maize molecular breeding with enhanced quality traits. 



Background 

The reversible protein phosphorylation is a fundamental 
mechanism that modulates many cellular functions includ- 
ing regulating developmental events and perceiving en- 
vironmental stimuli [1]. During phosphorylation, protein 
kinases (PKs) mainly phosphorylate serine (Ser), threonine 
(Thr) and tyrosine (Tyr) residues, while protein phospha- 
tases (PPs) can reverse this process by removing the phos- 
phate group. Based on substrate specificities, PPs can be 
categorized into two groups: Ser/Thr and Tyr phospha- 
tases. Recently, 1241 maize PK-encoding genes have been 
identified, with the data suggesting that maize PKs were 
implicated in diverse biological processes, such as develop- 
mental control and drought stress [2,3]. The development 
of a flowering plant both at the cellular and organismal 
level is a highly complex phenomenon. In spite of the 
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apparent autonomy, roots, leaves, and other structures 
perceive external signals and use different signaling mech- 
anisms to respond. Like functions of mammalian PP genes 
[4], plant PPs are expected to be the key components in 
signal transduction networks at distinct stages of plant de- 
velopment and in response to multiple abiotic signals. 
Specific inhibitor (okadaic acid and calyculin A) of Ser/ 
Thr phosphatases PPl and 2A-type protein phosphatase 
(PP2A) arrest root hair growth in a tiny range, and se- 
verely affect the shape of cells within the elongation zone 
and inhibit root growth. In addition, PP2A was associated 
with controlling microtubules length [5]. As reported in 
Arabidopsis, PLL4 and PLL5 (POL-like gene) regulate leaf 
development, but having no detectable functions within 
the meristem [6]. Flowering time is a major adaptive trait 
and an important selection criterion in plant breeding. 
Kim et al. [7] showed that a self-regulatory phytochrome 
kinase-phosphatase coupling is a key signaling component 
in the photoperiodic control of flowering in Arabidopsis. 
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Plants as sessile organisms are constantly challenged by 
a wide range of abiotic stresses such as salinity, drought, 
and cold. However, stresses are not necessarily a problem 
for plants because they have used several strategies to 
avoid or reduce the possible damage. Drought limits plant 
growth and result in a drastic decline in the photosyn- 
thetic yield due to osmotic stress-imposed constraint. In 
order to cope with drought stress, plants utilize various 
mechanisms, known as drought escape, drought avoidance 
and drought tolerance [8]. Previous study showed that 
Group A PP2Cs were negative regulators in ABA signaling 
pathway and acted as key regulators of desiccation toler- 
ance in land plants [9]. Salinity is one of the most severe 
environmental factors that greatly impacts plant develop- 
ment and restricts crop production. Furthermore, salt 
stress may cause water deficit, ion toxicity, nutrient imbal- 
ance, and oxidative stress. Interestingly, salinity tolerance 
is more likely to be a complex multigenic trait involving 
responses to cellular osmotic and ionic stresses. In Arabi- 
dopsis, over-expression of AtPP2CGl increased salt toler- 
ance, whereas its loss of function impaired salt tolerance 
[10]. Additionally, cold stress is another key environmental 
factor that limits the geographic distribution of plants. 
Cold acclimation has been described as the process in 
which plants adjust their metabolism during cold treat- 
ment. In spite of the fact that maize is generally sensitive 
to low temperatures [11], the extent of cold sensitivity is 
different within the maize germplasm. The cold signal is 
initially perceived by plasma membrane with the help of 
specific Ca^* channel proteins, membrane histidine ki- 
nases and some unknown sensors, which then activate the 
sophisticated cold-responsive signaling pathways in con- 
cert with plant hormone signaling, the circadian clock, and 
the developmental transition to flowering. Two PP2Cs, 
ABIl and AtPP2CA, which were reported as negative reg- 
ulators in ABA signaling, exhibited differences in their 
tissue-specific expressions as well as in temporal induction 
in response to chilling [12]. To improve the crop yields 
under abiotic stress conditions, it is crucial to understand 
the fundamental molecular mechanisms behind stress tol- 
erance in plants. 

One hundred twelve and one hundred thirty-two can- 
didate PP genes were identified in Arabidopsis and rice, 
respectively [13,14]. In Arabidopsis, Group A PP2Cs, 
ABIl and ABI2, are negative regulators of ABA signal- 
ing, whereas Group B are characterized to regulate 
mitogen-activated protein kinase (MAPK) signaling cas- 
cades. POL-type phosphatases, members of Group C 
PP2C, are involved in flower development. Independently, 
another PP2C member, kinase-associated protein phos- 
phatase (KAPP) is a singleton that regulates receptor-like 
kinases [15]. PTPs are critical partners for tyrosine- 
specific kinases in regulating the tyrosine phosphorylation 
status of many proteins. Tyrosine phosphorylation plays a 



role in MAPK signaling cascades. In particular, the level of 
tyrosine phosphorylation is determined by the balanced 
activity of protein tyrosine kinases (PTKs) and PTPs. 
Moreover, previous study showed that tyrosine phosphor- 
ylation is involved in phytohormone responses [16]. Com- 
pared with plants, other genomes which lack tyrosine 
kinases also have very few PTPs. Hence, identification of 
PTPs provides a stepping stone to better understanding of 
the functional significance of tyrosine phosphorylation in 
higher plants. 

Although recent advances in higher plants, only a small 
portion of the members in this family have been found to 
be components of signaling pathway in maize. This in- 
clude, ZmPP2Ca (ZmPP8), which is involved in the signal 
transduction in regulating response to drought stress [17]; 
ZmPP2C2 (ZmPP159), which plays a positive role in to- 
bacco cold resistance [18]; ZmPP2C (ZmPP76), which is 
implicated in Arabidopsis stress signal transduction [19]; 
ZMPP2 (ZmPP138), which was selected as a candidate for 
the catalytic subunit of phospho-pyruvate dehydrogenase 
phosphatase (PDP) [20]; ZmRIPl (ZmPPl), which acts as 
a chloroplast-to-nucleus signaling messenger, and is con- 
firmed to function in maize redox signaling [21]; ZmPPl 
(ZmPP65), whose function and regulation might be very 
similar to that of mammalian PPs [22]. 

In the post-genome era, the whole-genome sequen- 
cing, together with the global transcriptome profiling, 
such as microarray and RNA-seq, offer the opportunity 
to identify diverse gene families and to unravel the func- 
tions of genes involved in processes such as develop- 
mental regulation, disease resistance and abiotic stress. 
Maize is not only one of the most important food crops 
of the world, but also a model plant for study of the gen- 
etics in monocotyledons. In this study, we identified the 
full complement of PP genes in maize genome for evolu- 
tionary and functional analyses, which would help crop 
breeders to develop improved varieties. 

Methods 

Identification and characterization of ZmPPs 

The maize protein sequences were downloaded from the 
Maize Genome Sequence Project (http://ftp.maizesequence. 
org/release-5b/filtered-set/). To uncover all the members 
of PPs in maize, the predicted phosphatase sequences 
from Arabidopsis thaliana and Oryza sativa were used 
as query sequences to search against maize protein 
database using HMMER 3.0 software [23]. Firstly, we 
used hmmbuild tool from HMMER to build the PP hid- 
den Markov model (HMM) profile from the alignments 
of the known PP sequences. Secondly, hmmsearch pro- 
gram was separately applied to a search of all the protein 
sequences in maize protein database with the Pfam PP 
profiles (PF00481: PP2C, PF01451: LMWPc, PF00782: 
DSPc, PF00102: Y phosphatase) and PP HMM profile. 
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and the E-value cut-off was set to 1. Subsequently, each 
protein sequence was subjected to SMART (http://smart. 
embl-heidelberg.de/) and Pfam (http://pfam.sanger.ac.uk/) 
databases to ensure the presence of the catalytic domain. 
Proteins without a phosphatase catalytic domain were re- 
moved from the dataset. After eliminating the incomplete 
sequences manually, 159 sequences remained were finally 
identified as ZmPPs and renamed based on their loci on 
chromosomes. The subcellular localization of ZmPPs was 
predicted by WoLF PSORT (http://wolfpsort.org/). In 
addition, the molecular weights (MWs) and isoelectric 
points (pis) of ZmPPs were predicted by ExPASy Server 
tool (http://web.expasy.org/compute_pi/). 

Phylogenetic analysis and classification of the ZmPP gene 
family 

Catalytic domains of ZmPPs were used for multiple align- 
ments with the aid of MEGA5 software (http://www. 
megasoftware.net/) by employing ClustalW as the algo- 
rithm. Phylogenetic trees were constructed by neighbor- 
joining (NJ) algorithm of MEGA5. Bootstrap values from 
100 replicates are indicated at each node. Our PP classes 
were defined based on the the classification criterion sug- 
gested by Singh et al. [13,14]. 

Gene structure and intron/exon configuration 

To investigate the intron/exon structures of PP genes 
from maize and rice, we collected useful information 
from genome annotations of maize and rice from the 
Maize Genome Sequence Project and TIGR database, re- 
spectively. Both DNA sequences and the corresponding 
coding sequences were loaded into the Gene Structure 
Display Server (http://gsds.cbi.pku.edu.cn/). For a better 
visualization and comparison, the 5' untranslated region 
(UTR) sequences were removed beforehand. 

Chromosomal localization, gene duplication and synteny 
analysis 

To determine the location of ZmPP genes on 10 chro- 
mosomes, we extracted data concerning gene positions 
from the Maize Genome Sequence Project. MapDraw 
2.2 was used to visualize the locations of ZmPP genes on 
maize chromosomes with physical distances in million 
bp (Mb). We adopted a group-specific color strategy to 
mark each gene for better visualization and recognition. 
Tandemly duplicated genes were defined as closely re- 
lated genes on a single chromosome, with no more than 
ten intervening genes [2]. In contrast, segmentally dupli- 
cated genes were detected by CoGe SynMap program 
(http://genomevolution.org/CoGe/SynMap.pl). In order 
to determine the syntenic regions between the rice and 
maize genomes, the positions of syntenic regions from 
these two genomes were collected from the results cal- 
culated by Synteny Mapping and Analysis Program 



(SyMAP 4.0) [24]. The segmental duplications and syn- 
tenic regions were finally visualized using Circos 0.62 
(http://circos.ca). 

The number of nonsynonymous substitutions per non- 
synonymous site (Ka) and the number of synonymous 
substitution per synonymous site (Ks) of duplicated 
genes were calculated by DnaSP 5.0 [25]. The ratio of 
non-synonymous to synonymous nucleotide substitu- 
tions (Ka/Ks) between paralogs was analyzed to explore 
the mechanism of gene divergence after duplication. The 
dates of the duplication events were calculated by the 
equation T= Ks/2A x 10"^ Mya, the A = 6.5 x lO"** [26]. 

Expression analyses by microarray 

To analyze the expression patterns of ZmPPs, the tran- 
scriptome data of genome-wide gene expression atlas of 
maize inbred line B73 were downloaded from the database 
PLEXdb (http://www.plexdb.org/) with the accession num- 
ber ZM37. For drought treatments, microarray data were 
obtained from NCBI Gene Expression Omnibus (GEO) 
with accession number GSE16567. For each treatment, 
signal intensity values were normalized and biological rep- 
licate samples were averaged to generate log2 expression 
values for each gene. Log2-transformed values were loaded 
into R (15.2) and Bioconductor for expression analysis 
(http://www.bioconductor.org/). Limma package was ap- 
plied to data processing, and heatmaps representing log2- 
transformed values were generated with the gplots package 
[27]. A hierarchical clustering algorithm was applied to de- 
termine similar pattern in expression profiles. Microarray 
data of cold stress were extracted from the ArrayExpress 
database under the experiment accession number E- 
MTAB-1315. For drought, differentially expressed genes 
were selected under a very stringent cutoff, with 0.05 P- 
value and fold change value of 2, whereas the response 
must be at least 2.8-fold in cold stress. 

Differential gene expression profiles based on RNA-seq 

To generate the expression profiles of ZmPPs among dif- 
ferent organs and development stages, the RNA-seq data 
were downloaded from the NCBI Short Read Archive with 
accession number SRP010680 (http://www.ncbi.nlm.nih. 
gov/sra/). This RNA-seq data analyzed 18 selected tissues 
representing five organs from maize inbred line B73. Nor- 
malized gene expression values were estimated by frag- 
ments per kilobase pair of exon model per million 
fragments mapped (FPKM). Finally, log2-transformed 
FPKM values from 18 tissues were used to draw heatmaps. 
To avoid taking the log of a number less than 1, all such 
FPKM values were replaced by 1. Similar to microarray 
analysis, both limma and gplots packages were used to 
generate the heatmaps. 

RNA-seq data of three types of roots (CR, PR, and SR) 
under salt stress were fetched from GEO database at 
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NCBI under the corresponding accession number 
GSE53995. Normalized signal intensity values were log- 
transformed for further analysis. Moreover, 2-fold 
change was used to test the differentially expressed 
genes. The up- or down-regulated genes in any tissues 
were calculated from the average of log-transformed 
normalized signal values. 

Plant materials and stress treatments 

Plants of maize were grown in experimental plots with a 
1:1:1 mix of peat:vermiculite:perlite. A growth room with 
controlled environmental conditions (15 h light/25°C, 

9 h dark/20°C, relative humidity of 55%) was used for 
growing plants. The seedlings were watered daily for 
about 2 weeks. For drought treatment, 14-day-old seed- 
lings of maize inbred line Ye478 and Han21 were ex- 
posed to drought stress for 3 or 4 days, while controls 
were well watered throughout this period. At the end of 
treatment period, the volumetric water content of the 
soil (measured at the depth of approximately 5 cm from 
the top of soil level) was approximately 12.5% for the 
drought-stressed plants, compared with approximately 
30% for the control plants. For salt treatment, maize in- 
bred line B73 was selected. Seedlings were watered with 
200 mM NaCl for 24 h, after which all stress-treated and 
control seedlings were harvested, separated into shoots 
and roots and stored at -80°C. Each treatment was repli- 
cated three times with 5 to 7 seedlings. For cold treatment, 
kernels of maize inbred line Huangzao4 and Moll7 were 
germinated for 3 d in darkness at 25°C. Seedlings 
were transferred to pots containing Knop's solution and 
were grown in a growth chamber (photoperiod: 14/ 

10 h day/night, light irradiance: 250 |imol quanta m s 
temperature: 24/25°C and relative humidity 60/80%). 
When the third leaf was fully developed, half of the plants 
were transferred to 8/6°C (day/night) without changing 
the other conditions; the other half were grown in the 
same conditions as before (control plants). 

Seeds of Arabidopsis thaliana were sterilized with 2% 
bleach for 20 min, plated on Murashige and Skoog 
medium (MS; Sigma-Aldrich), chilled at 4°C for 3 d, and 
transferred to a growth room at 19 to 23°C with a 16 h/ 
8 h light/dark photoperiod. 

RNA extraction and real-time PCR 

To validate the RNA-seq data for selected genes under 
abiotic stresses, real-time PCR was carried out. The 
leaves and roots of the maize seedlings were gathered, 
immediately frozen in liquid nitrogen and then stored at 
a temperature of -80°C. Subsequently, total RNA was 
isolated from selected tissues using Trizol reagent, and 
the use of DNase I treatment to remove any genomic 
DNA contamination from RNA samples. About 2 |ig of 
total RNA was reverse transcribed using a Takara RNA 



Kit to generate the single-stranded cDNA. Real-time PCR 
was performed on iCycler iQ5 Multicolor real-time PCR 
detection system (Bio-Rad) by using the Power SYBR 
Green PCR Master Mix (AP Plied Biosystems). Two bio- 
logical replicates of each sample were used and three 
technical replicates were performed for each biological 
replicate. The thermal cycling conditions are as follows: 
50°C for 2 min, 95°C for 10 min, 40 cycles of 95°C for 
15 s, and 60°C for 1 min. The relative quantification 
method 2' ' was used to evaluate quantitative variation 
between replicates. For drought and cold stresses, maize 
ubiquilin-1 (UBQl) gene was used as internal control, 
while ACTIN7 was an internal control for salt stress. 

For fptpl mutant analysis, total RNA was extracted 
from rosette leaves of 40 day-old plants. ACTIN8 was 
used as an internal control. The real-time PCR experi- 
ment was repeated for three independent biological rep- 
licates. Finally, the ABi7500 real-time PCR system with 
the RealMasterMix (SYBR Green I) (Takara) was used 
for performing real-time PCR. Primer pairs used in this 
study were listed in Additional file 1. The fptpl mutant 
was obtained from Institute National de la Recherche 
Agronomique (INRA). 

Over-expression of FPTPl In Arabidopsis and analysis of 
ABA synthetase gene expression 

To generate Arabidopsis FPTPl over-expression trans- 
genic plants, the coding sequence of FPTPl were amplified 
and inserted into the modified pCAMBIA1300-Super vec- 
tor under control of 35S promoter. These vectors were 
transformed into Agrobacterium strain GV3101 and subse- 
quently introduced into the Arabidopsis wUe-type plants. 

With the aim of identifying the contribution of PTPs in 
ABA accumulation, RNA was extracted using an RNeasy 
Plant Minikit (QIAGEN), and cDNA was synthesized 
using M-MLV Reverse Transcriptase (M1701) (Promega). 
1 |ig RNA was mixed thoroughly with 1 |ig Oligo (dT)i5, 
and treated at 70°C for 5 min. Subsequently, the following 
were added: 5 M-MLV buffer, 1.25 |iL dNTP Mbc, 
1 i^L M-MLV, and 0.6 ^iL RNase inhibitor. The resulting 
25 i^L reaction was incubated for 1 h at 37°C. PCR primers 
were showed in Additional file 1. PCR cycles included 
94°C 5 min, (94°C 45 s, 58''C 45 s, 72°C 1 min) 35 cycles, 
72°C 7 min. 

Subcellular location of ZmPP1 and estrogen treatment of 
ZmPPl 

For subcellular localization of ZmPPl, its cDNA se- 
quence was cloned into pEZS-NL transient expression 
vector to generate pEZS-NL-ZmPPl. Particle bombardment- 
mediated transient expression of ZmPPl-GFP fusion pro- 
tein in onion epidermal cells were investigated by laser 
scanning confocal microscopy. As for GUS staining, the 
cDNA of ZmPPl was cloned into pSuper-1391 vector and 
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then introduced into the GV3101. Floral dip method can 
efficiently generate the Arabidopsis transgenic plants. 
Then the transgenic plants were screened on solid plates 
containing 50 mg/mL hygromycin. 

We have developed an estrogen receptor-based chemical- 
inducible system for ZmPPl transgenic Arabidopsis. The 
amplified ZmPPl cDNA sequence was cloned into pER8, 
of which the pER8 vector was provided by N.-H.Chua and 
B. Ulker [28]. Transformations of Arabidopsis were per- 
formed by the floral dip method using GV3101. Over- 
expression of ZmPPl in the transgenic lines was induced 
by spraying the seedling with 5 or 10 i^M |3-estradiol. 

Promoter activity and computational prediction of miRNA 
protein-coding target transcripts 

To find out putative cis-acting elements in the promoters 
of putative ZmPP genes, the 1,000 bp upstream promoter 
regions of all ZmPPs were used to search for known 
stress-responsive ds-elements by PlantCARE database 
(http://bioinformatics.psb.ugent.be/webtools/plantcare/html/). 
Target-align [29], a miRNA target prediction tool, was 
used to predict putative miRNA target genes. For the 
maize protein-coding transcripts, the predicted cDNA of 
the longest consensus maize transcript were used. A total 
of 321 mature miRNA sequences downloaded from miR- 
Base database (release 20) [30] were reverse complemen- 
ted and matched against the indexed maize transcript 
database. With the default setting, a score cut-off of > 75 
and mismatch < 4 nucleotides were applied in the 
prediction. 

Results and discussion 

Identification of PPs in maize genome 

Previously, 112 and 132 PP genes were identified from 
Arabidopsis and rice, respectively [13,14]. To explore the 
occurrence and size of the PP family in maize, HMMER 
program was used to search against the maize proteo- 
mics database. Further, the presence of catalytic domain 
was confirmed by SMART and Pfam databases. From 
maize genome, we identified a total of 159 putative PP- 
encoding genes following the removal of those se- 
quences with an incomplete catalytic domain, suggesting 
that the size of maize PP family was larger when com- 
pared with that in Arabidopsis and rice. For simplicity, 
the 159 ZmPP genes were renamed from ZmPPl to 
ZmPP159 according to their exact positions on maize 
chromosomes 1-10 and from top to bottom. The identi- 
fied PP genes in maize encode proteins ranging from 
131 to 1264 amino acids (aa). With the WoLF PSORT 
analysis, most PPs were predicted to locate in nucleus 
and chloroplast as well as other organelles. ExPaSy ana- 
lysis suggested that the ZmPP proteins had large varia- 
tions in isoelectric point (pi) values (ranging from 4.46 
to 9.8) and molecular weights (ranging from 14,245.15 



to 90,157.14 Da). All of the related information on 
ZmPPs is listed in the Additional file 2. 

Phylogenetic analysis and classification of the ZmPPs 

The catalytic domains of PPs are diverse even if the cata- 
lytic cores of all eukaryotic PKs share extensive similar- 
ities in both primary and three-dimensional structures. 
With the aim of analyzing the phylogenetic relationships 
of 159 ZmPPs among maize, phylogenetic tree was car- 
ried out based on the catalytic domains by NJ algorithm 
(Figure 1). Based on the sequence similarities of the 
ZmPP gene sequences and the classification of PP genes 
from Arabidopsis and rice, the phylogenetic tree divided 
the maize PP family proteins into three major classes, 
namely, PP2C, PTP and PPP (PP2A). PP2C was evolu- 
tionarily related to the major class of Ser/Thr PPs with 
104 genes, whereas PTP and PP2A comprised of 29 and 
26 members, respectively. Interestingly, all the members 
of PP2C class formed a single major group and this 
group could be further classified into 11 subgroups ex- 
cept some ambiguous branches. Each subgroup repre- 
sents a subfamily of PP2C and is designated from A-K 
according to Xue et al. [31] and Singh et al. [14]. Group 
J PP2C did not include any rice and maize PPs but con- 
tained only members from Arabidopsis, suggesting that 
they might have been lost both in rice and maize after 
being originated from common ancestor. 

In terms of the previous studies [1,31], plant lacks typ- 
ical tyrosine kinase, and most plant PPs (67%) have Ser/ 
Thr specificity. It seems that plants prefer to use Ser/Thr 
as substrates for PKs and PPs, in sharp contrast with the 
case in human [32]. Three genes belonging to LMWP 
(ZmPP14), PTPL (ZmPP47) and CDC25 (ZmPP145) clas- 
ses were positioned separately. PP2As, which were ubiqui- 
tous enzymes in eukaryotes, formed another single major 
class. Phylogenetic tree obtained from Arabidopsis, rice 
and maize showed very similar topologies and subfamily 
organization with individual maize tree (Additional file 3). 
It is noteworthy that the number of maize PP genes was 
generally overrepresented than that of Arabidopsis and 
rice in almost all classes. 

Gene structure analysis and intron/exon arrangement 

As a evolutionary relic, it is noteworthy that intron/exon 
arrangement carries the imprint of the evolution of a gene 
family, and may provide insight into their evolutionary 
mechanisms underlying the origin of gene families [33]. 
With the purpose of understanding the structure diversity 
and evolution of ZmPP gene family, we analyzed the in- 
tron/exon structure, intron position and phase. A detail il- 
lustration of the intron/exon structure was shown in 
Additional file 4. It was found that the great majority of in- 
trons were positioned in the coding sequences of ZmPP 
genes. Notably, in PP2C class, the intron/exon structures 
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Figure 1 Phylogenetic relationships of ZmPPs. An un-rooted NJ tree is made based on the catalytic domain sequences of ZmPPs. The whole 
PP gene family is divided into different classes, PP2C, PP2A, PTPs, PTPL, CDC25, and LMWP. PP2C class is further subdivided into different classes 
(A-K). Scale bar represents 0.1 amino acid substitutions per site. 



were different between subclasses and correlated with the 
architecture of protein sequence and phylogenetic analysis. 
However, it was notable that a similar gene structure 
was found in each group, while 4 members of the 104 
ZmPP2C genes {ZmPP39, -112, -91, and -142) had no in- 
trons. This indicated that the intron patterns, which 



correlate well with the phylogentic clades, strongly support 
their close evolutionary relationships among the ZmPP2C 
genes within the same subclass. Among those having in- 
trons, the number of introns within the open reading 
frame (ORF) ranged from 1 to 12, showing a great differ- 
ence in the ZmPP2C class. A great degree of variation in 
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the number of introns exists in groups E, G and K, while 
the number of introns in the rest groups changes within a 
small range, mostly from 1 to 3. In rice, similar results of 
gene structure analysis were given in Additional file 5, in- 
dicating that similar intron loss or gain events occurred 
during expansion and structural evolution of 2C-type pro- 
tein phosphatase may exist in monocotyledon. 

In PTPs, introns are not equally distributed across 
family. Additionally, we found that nearly all members 
(expect ZmPP73) have introns, while the number of in- 
trons varies widely within PP2A group. We further ana- 
lyzed the exon/intron configuration of the paralogous 
pairs in PP gene family to obtain traceable intron gain or 
loss information (Additional file 5). Despite the struc- 
tural conservation of intron/exon arrangement found in 
some paralogous pairs, others exhibited extensive vari- 
ation. What's more, growing evidence has showed a 
functional link between the structural diversity of gene 
members and the evolution of multiple gene families, 
while intron loss or gain can be an important step in 
generating structural diversity and complexity [34]. 

Chromosomal location and gene duplication of ZmPP 
genes 

To determine the genomic distribution of the ZmPP genes 
and study their evolution in the context of whole genome 
duplication, we mapped each gene on maize chromo- 
somes based on their corresponding coordinates. As 
shown in Figure 2, ZmPPs are almost unevenly distributed 
across the ten chromosomes. All members of class PP2C 
are distributed across 10 maize chromosomes. In contrast, 
no PTPs and PP2As were found to be located on chromo- 
some 7 and 10, respectively (Figure 2). 

Gene duplication is prominent in eukaryotic evolution, 
because duplicated genes provide the raw materials for 
evolving a new gene, which in turn facilitate the gener- 
ation of novel gene functions. Compared with most other 
eukaryotes, plants have substantially higher gene duplica- 
tion rates in the adaptive response to environmental stim- 
uli [35]. Previous study revealed that maize genome 
underwent a whole-genome duplication event (WGD) that 
resulted in the presence of two sub-genomes and -30% of 
genes retained copies in both sub-genomes [36]. The 
duplication mechanism including transposition, seg- 
mental duplication, replicative transposition, or even 
WGD underlying gene family expansion are similar be- 
tween eukaryotes. To clarify the probable relationship be- 
tween ZmPP genes and potential gene duplication within 
the genome, we analyzed the occurrence of tandem dupli- 
cation and large-scale segmental duplication during the 
evolution. In all, 15 (-9.4%) ZmPP genes were found to be 
tandem repeats with a maximum number of ten interven- 
ing spacer genes (Table 1). These 15 ZmPP genes were 
represented in 9 distinct tandem duplicated gene clusters. 



with one cluster containing 3 tandem genes and the rest 
of clusters possessing 2 tandem genes. It is notable that 
genes in one duplicated cluster may come from different 
group. For instance, ZmPP144 was grouped into PP2C, 
but ZmPP145 was a member of PTPs. Hence, it is specu- 
lated that tandem duplications have an indispensable role 
in the evolution of maize PP gene family. Additionally, 
previous study showed that most plants were diploidized 
polyploids and retained numerous duplicated chromo- 
somal blocks in their genomes, segmental duplication oc- 
curred most frequentiy in plants [37]. To evaluate the 
contribution of large segmental duplications in the expan- 
sion of ZmPP-encoding gene family, CoGe SynMap pro- 
gram was used to detect segmental duplicated genes. We 
noted that 38 segmental duplication events in maize were 
observed, which were originated from a polyploidy event 
occurred around 15 to 20 million years ago (Table 1). 
Interestingly, all genes in segmental duplicated pairs were 
from the same classes. Thus, we speculated that segmental 
duplications exclusively contributed to the expansion of 
PP gene family in the maize genome. Taken together, our 
analyses suggested that large-scale segmental duplication 
appear to have prominently contributed to the current 
complexity of the maize PP family. It is consistent with 
the duplicative mechanism of PP gene family in rice [14]. 

To illuminate the divergence after gene duplication, the 
ratios of Ka/Ks were estimated for all 48 duplicated pairs. 
Usually, Ka/Ks > 1 means positive selection; Ka/Ks = 1 
means neutral evolution; while Ka/Ks < 1 means purifying 
selection [38,39]. Ka/Ks ratios of most duplicated pairs, re- 
gardless of whether they are orthologs and paralogs, were 
less than 1, subjecting to purifying selection. Nevertheless, 
there are 5 segmentally duplicated pairs suffered positive 
selection with Ka/Ks ratios greater than 1. Among these 
duplicated pairs, the average Ka/Ks value of segmental du- 
plication was very similar to those in tandem duplication 
(-0.4), indicating that most duplicated ZmPP genes are 
under strong purifying selection pressure. The duplication 
events for the segmental duplicated genes were estimated 
to occur approximately 24 Mya, while that in tandem du- 
plicated pairs were estimated to occur approximately 37 
Mya, indicating that the segmental duplication events 
might occur before the emergence of tandem duplication 
events. Therefore, the segmental and tandem duplication 
events played essential roles in evolutionary expansion of 
PP family in maize. 

Expression profiles of maize PP genes and their potential 
functions in different tissues and developmental stages 

Recent advances in functional analyses of PPs have re- 
vealed the importance of some PP genes in the life cycles 
of both Arabidopsis and rice [14,40], however, their roles 
in maize remained unclear. Aim for achieving gene ex- 
pression patterns of ZmPPs in diverse growth phases, we 
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Figure 2 Chromosomal localization of ZmPP genes on 10 chromosomes of maize. Respective chromosome numbers are written at the top. 
Genes belonging to three classes have been marl<ed by different colors (PTPs means protein tyrosine phosphatases). Dashed lines join the genes, 
lying on duplicated segments of the genome. Tandemly duplicated genes are joined with vertical lines. 



downloaded the genome-wide expression data of maize 
inbred B73 from PLEXdb and analyzed the expression 
profiles of ZmPPs across the 60 different developmental 
stages of 11 organs. Following whole-chip data process- 
ing, the log2 values of 154 PP genes were showed as 
heatmap (Additional file 6). From Additional file 6, we 



found that ZmPPs were implicated in a wide range of 
plant development processes, including root develop- 
ment, leaf development gradient and flower induction. 
In addition to conducting microarray analysis, RNA-seq 
data with 18 tissues representing distinct stages of maize 
plant development were selected for analysis. A total of 
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Table 1 The Ka/Ks ratios and estimated absolute dates for the duplication events between the duplicated ZmPP genes 



Duplicated pair 


Ka/Ks 


Date (Mya) 


Duplicate type Purifying selection 


Group 


ZmPP53/ZmPP132 


0.2091 


1 3.02307692 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP56/ZmPP130 


0.35128 


13.86153846 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP107/ZmPP127 


0.26416 


18.2 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP134/ZmPP66 


0.10751 


86.57692308 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP8/ZmPP147 


0.13034 


20.53846154 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP7/ZmPP148 


0.39908 


10.06153846 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP24/ZmPP77 


1 .29469 


5.507692308 


Segmenta 


No 


PP2C/PP2C 


ZmPP30/ZmPP155 


0.15067 


12.1 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP4/ZmPP151 


0.10791 


16.53846154 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP21/ZmPP80 


0.14157 


21.40769231 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP28/ZmPP92 


0.16562 


67.48461538 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP83/ZmPP99 


0.08266 


14.23846154 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP41/ZmPP116 


0.11765 


16.73846154 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP27/ZmPP152 


1 .60306 


16.60769231 


Segmenta 


No 


PP2C/PP2C 


ZmPP152/ZmPP158 


2.32622 


11.17692308 


Segmenta 


No 


PP2C/PP2C 


ZmPP158/ZmPP27 


0.21536 


1 232307692 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP72/ZmPP86 


0.11398 


16.06153846 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP76/ZmPP84 


0.09925 


15.34615385 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP122/ZmPP129 


0.19247 


65.58461538 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP95/ZmPP65 


0.07717 


30.7 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP9/ZmPP141 


0.18169 


11.17692308 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP29/ZmPP156 


0.16686 


13.23076923 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP42/ZmPP118 


0.10881 


14.13846154 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP85/ZmPP100 


1 .28064 


17.76153846 


Segmenta 


No 


PP2C/PP2C 


ZmPP23/ZmPP78 


0.17848 


23.23076923 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP37/ZmPP110 


0.25391 


1 5.238461 54 


Segmenta 


Yes 


PP2C/PP2C 


ZmPP33/ZmPP68 


0.18219 


16.76153846 


Segmenta 


Yes 


DSP/DSP 


ZmPP52/ZmPP131 


0.12905 


14.00769231 


Segmenta 


Yes 


DSP/DSP 


ZmPP54/ZmPP106 


0.24256 


70.78461538 


Segmenta 


Yes 


DSP/DSP 


ZmPP67/ZmPP133 


0.33671 


44.04615385 


Segmenta 


Yes 


DSP/DSP 


ZmPP94/ZmPP133 


2.1275 


11.94615385 


Segmenta 


No 


DSP/DSP 


ZmPP94/ZmPP101 


1.8272 


13.8 


Segmenta 


No 


DSP/DSP 


ZmPP101/ZmPP133 


0.13011 


13.36153846 


Segmenta 


Yes 


DSP/DSP 


ZmPP11/ZmPP50 


0.05879 


15.17692308 


Segmenta 


Yes 


PP2A/PP2A 


ZmPP11/ZmPP81 


0.01588 


105.6230769 


Segmenta 


Yes 


PP2A/PP2A 


ZmPP19/ZmPP81 


0.05683 


14.75384615 


Segmenta 


Yes 


PP2A/PP2A 


z_iIlrrZZ/Z_lllrr/y 


n nn7Q9 
u.uu/ yz 


1 ^01^07 f^Q 

\ j.jyzju/ oy 


Segmenta 


Yes 


PP9A/PP9A 


ZmPP75/ZmPP88 


0.05525 


9.884615385 


Segmenta 


Yes 


PP2A/PP2A 


ZmPP2/ZmPP3 


0.10884 


8.269230769 


Tandem 


Yes 


PP2C/PP2C 


ZmPP2/ZmPP4 


0.12221 


5.853846154 


Tandem 


Yes 


PP2C/PP2C 


ZmPP3/ZmPP4 


0.1199 


6.030769231 


Tandem 


Yes 


PP2C/PP2C 


ZmPP5/ZmPP6 


0.96739 


84.93076923 


Tandem 


Yes 


PP2A/PP2C 


ZmPP60/ZmPP61 


0.36275 


41.66923077 


Tandem 


Yes 


PP2C/PP2C 


ZmPP119/ZmPP120 


0.77549 


49.03076923 


Tandem 


Yes 


PP2A/DSP 
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Table 1 The Ka/Ks ratios and estimated absolute dates for the duplication events between the duplicated ZmPP genes 

(Continued) 

ZmPP138/ZmPP139 0.43648 2361538462 Tandem Yes PP2C/PP2C 

ZmPP141/ZmPP142 0.57337 85.86153846 Tandem Yes PP2C/PP2C 

ZmPP144/ZmPP145 0.68466 50.03076923 Tandem Yes PP2C/DSP 

*Ka/Ks < 1 means purifying selection. 



152 ZmPPs were reliably detected, while the remaining 7 
members were undetected (Figure 3A, Additional file 7). 
We also calculated the CV value (CV = S/Xmean> where S 
represents the standard deviation and X^ean indicates 
the mean expression of a gene across all the tissues) of 
each gene to elucidate the expression profiles of differ- 
entially expressed genes [41]. Based on their expression 
levels and CV values, these detected PP genes were cate- 
gorized into two lineages. Lineage I showed higher ex- 
pression levels, but lower CV values compared with 
lineage II. As shown in Additional file 7, genes with 
higher CV values displayed higher fluctuation in their 
expression levels. 

Recently, a systematic survey of morphology and hist- 
ology of young seedling root in maize hybrids and their 
parental inbred lines has been performed [42] . Before this, 
most of our knowledge about the complexity of auxin 
signaling in regard to root development was based on 
the molecular and physiological analyses in Arabidopsis, 
which laid the groundwork for the discovery of parallel 
pathways using maize as a model for monocot plants. To 
better understand the role of auxin in maize root develop- 
ment, we searched for some key factors in the signaling 
pathways (Additional file 8). In comparison with Arabi- 
dopsis and rice, less development-regulated genes and 
their functions in maize have been characterized in detail. 
In our study, identification of key genes involved in root 
development signaling pathway in maize, such as auxin 
transporters (PIN) and auxin response factors (ARF) were 
performed using HMM-based searches. All genes involved 
in this developmental signaling pathway were showed in 
Additional file 9. It has been reported that PP2As are re- 
quired for asymmetric auxin distribution in seedling roots 
of Arabidopsis [43]. During male and female inflorescence 
differentiation and kernel development, the transcripts of 
ZmPIN gene were showed to have overlapping expression 
domain in the root apex [44]. Moreover, the inhibition of 
PP2A leads to a preferential apical PIN targeting [43]. In 
our microarray data, ZmPP65 was relatively high expressed 
in root, indicating that ZmPP65 might play an important 
role in auxin-induced root development. Interestingly, an- 
other two genes {ZmPP67 and ZmPP94), which belong to 
PP2C class, exhibited high expression levels in root. Thus, 
we hypothesize that except for PP2A, other classes of PP 
genes may participate in root development as well, such as 
PP2C. Moreover, expression patterns of other PP genes 



involved in maize root developmental pathway were ana- 
lyzed as well (Figure 3A, Additional file 9). Most of them 
were highly expressed in root, which were in accordance 
with the results that other PP genes except PP2As play 
roles in root development as well. 

It has been reported that maize leaves are characterized 
by clear proximal/distal domains as seen in the morpho- 
logical differences of the sheath and blade. In our research, 
it was noteworthy that two genes (ZmPPMS and ZmPPSS), 
belonging to lineage I, were found to display high expres- 
sion levels at V5_Tip of stage-2 Leaf, while ZmPPSS was 
highly expressed at V9_Eighth leaf It is speculated that 
ZmPPMS and ZmPPSS may play critical roles in early leaf 
developmental stages (V5 and V9). To provide a broader 
insight into the maize leaf development, four representa- 
tive zones of the leaf blade were selected for deeper tran- 
scriptomic analysis (Additional file 10) [45]. Of 159 ZmPP 
genes, a total of 88 genes were expressed at all four sam- 
pled sections of the maize leaf Most of them were up- 
regulated in all four regions of maize leaf. Especially, 
ZmPP47 showed increased transcript abundance in the 
basal (base, 1 cm above the leaf three ligule) and transi- 
tional zones (-1 cm, 1 cm below the leaf two ligule), while 
ZmPPSS and ZmPP 133 were highly expressed in the ma- 
turing (+4 cm, 4 cm above the leaf two ligule) and mature 
zones (tip, +1 cm below the leaf three tip), which might 
participate in the transcriptional regulation of carbon fix- 
ation pathways. Arabidopsis PLL4 (a member of PP2C 
class) was reported to play a part in the regulation of leaf 
development, while its homolog in maize (ZmPP24) 
showed higher expression in the maturing and mature 
zones as well. Strikingly, 4 genes (ZmPP4S, -113, -120 
and -146) were obviously down-regulated in the maturing 
and mature zones with their functions in leaf development 
remain unknown. Overall, genes with specific expression 
patterns in this study might play critical roles in leaf 
development. 

Morphogenesis is referred to as a process in develop- 
mental biology in which organisms develop their charac- 
teristic shapes. Furthermore, flowering and architecture 
development are two key processes in plant morphogen- 
esis, which are of particular significances in crop. Modu- 
lation of the life cycle by photoperiodic signaling is a 
fundamental mechanism for plant to cope with the sea- 
sonal changes. As reported by Dong et al. [46], both en- 
vironmental and endogenous factors control the shoot 
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Figure 3 RNA-seq analysis of ZmPPs in 18 selected tissues and a genetic regulatory network for flowering time control in maize. 

A Expression profiles of 152 detected ZmPPs in 18 tissues from different developmental stages. B Real-time PCR analysis of CO and FT in fptpl 
mutant. C Plant images of fptpl transgenic plants. Left, wild type plant; right, FPTPl overexpression transgenic line. D Overview of flowering time 
pathway in maize. 
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apical meristem (SAM) transition from the vegetative to 
the reproductive stage. In Figure 3A, ZmPP67 was 
expressed at a relative high level in both root and SAM, 
suggesting that ZmPP67 might participate in floral tran- 
sition. To validate the exact role of ZmPP67 in maize 
flowering time pathway, the function of Arabidopsis 
FLOWERING ASSOCIATED PTPasel (FPTPl, Gene- 
bank accession number: FJ605097), which is homolo- 
gous to ZmPP67, was analyzed by mutant analysis. CO, 
a circadian-regulated gene, and FLOWERING LOCUS T 
{FT), belonging to the BBX and PEBP family, respectively, 
contribute to the photoperiod flowering control in plant 
[47]. As showed in Figure 3B, mutation in FPTPl signifi- 
cantly promoted FT expression but had no obvious effect 
on CO expression. Furthermore, over-expression of FPTPl 
caused delayed flowering (Figure 3C). Thus, we suspected 
that ZmPP67 may function in the upstream of FT. 

It was reported that many of the flowering time path- 
ways in Arabidopsis and rice are conserved in maize [47]. 
Through a genetic and molecular model for flower devel- 
opment in Arabidopsis, a possible flowering time pathway 
in maize was proposed (Figure 3D). Light is the most im- 
portant environmental cue implicated in the regulation of 
flowering time in plants. In maize, a reduced light re- 
sponse is associated with the development of early flower- 
ing inbred lines [48]. Phytochromes and cryptochromes 
are the primary red/far-red and blue light photoreceptors, 
respectively, with three pairs in our study, PhyBl/2, 
PhyAl/2, and CRYl/2. In our expression profiles of genes 
involved in flower induction (Figure 3D), PhyA, PhyB and 
CRYl/2 were expressed slightly higher levels in leaves 
compared with other tissues. Particularly, expressions of 
PhyBl and PhyB2 are highly similar to that of PhyAl and 
PhyA2 genes. Previous study showed that PhyB regulates 
CO protein, but not mRNA level [49]. CO coordinates 
light and clock input in leaves to trigger the expression of 
florigen gene FT. Consequently, we suspected that the 
effect of CO on FT expression must result from a post- 
transcriptional regulation. DLFl, homologous to Arabi- 
dopsis FD, form a complex with ZCN8 (maize FT protein) 
to activate downstream floral organ identity genes, such as 
ZMM5, the homolog of Arabidopsis SOCl in maize. Then, 
APl/3, PI and AG were the downstream of LEAFY in 
Arabidopsis, and the homologs in maize were identified in 
this study as well as their expression patterns were ana- 
lyzed. It was showed that DLFl, ZMMS, ZFL2 (LEAFY), 
GL15 {AP2), Silkyl {AP3), and ZMM16 (PI) were relatively 
highly expressed in SAM compared with other tissues, 
while other genes in the downstream of LEAFY, including 
ZAPl (API) and ZAGl/2 (AG) had lower expression. All 
genes involved in the flowering time signaling pathway 
were showed in Additional file 11. 

In addition to RNA-seq analysis in 18 representative 
tissues during maize development, another RNA-seq 



study of 12 maize reproductive tissues that represent 
male, female and developing seed was carried out. 
Overall, -96% (153) of ZmPPs were detected in reproduct- 
ive transcriptome (Additional file 12). It was noteworthy 
that ZmPP37, -76, -81, -84, -90, and -110 were signifi- 
cantly highly expressed in pollen, notably, ZmPP76 and 
ZmPP84 were specially expressed in anther. Additionally, 
ZmPPlS, -63 and -97 had relatively high expression 
values in both anther and pollen. Hence, genes with ex- 
pression restricted to specific tissues within this study are 
strong candidates for further analysis to unravel the mech- 
anisms of vegetative and reproductive growth. 

Expression profiles of ZmPPs under abiotic stresses 
Expression analysis of ZmPP genes under salt stress 

Plants encounter a wide range of environmental insults, 
which pose serious threat to crop production. Among all 
the abiotic stresses, salinity represents the major constraint 
for agricultural productivity. Maize is a salt-sensitive crop 
plant and responds to salt stress via complex changes at 
both the transcriptional and post-transcriptional levels. 
Plant roots are highly sensitive organ and induce primary 
response to salt stress. To better understand the different 
responding mechanisms to salt stress in maize root, exam- 
ination of three types of roots (CR, PR, and SR) under salt 
treatment was performed. As shown in the Figure 4A, 
three specific genes were significantiy and differentially 
expressed under salt stress, namely, ZmPP149, ZmPP66, 
and ZmPP127. Surprisingly, all of them are the members 
of PP2C class, and ZmPP66 is a member of Group G 
PP2C. It has been demonstrated that AtPP2CGl (Arabi- 
dopsis thaliana protein phosphatase 2C G Group 1) 
can positively regulate salt-tolerance of Arabidopsis in 
ABA-dependent manner [10]. Intriguingly, ZmPPlO, which 
is homologues to AtPP2CGl, showed significantly up- 
regulated expression in PR. In addition to these three spe- 
cially expressed genes (ZmPP149, ZmPP66, and ZmPP127), 
ZmPP107, a member of Group A PP2C, was significandy 
up-regulated in CR and SR, indicating that Group A PP2C 
may play a part in salt stress except for its role in drought 
stress condition. 

To clearly understand the roles of ZmPPs in salt stress- 
induced signaling pathway, salt stress signaling and mech- 
anisms of maize salt tolerance were analyzed. Salt stress 
involves both osmotic stress and ionic stress that limit the 
efficiency of crops and especially quantity and quality of 
their metabolic (secondary plant products) products. Os- 
motic stress can rapidly increase ABA biosynthesis, and 
regulate ABA-dependent stress response pathway [50]. 
The expression patterns of key genes involved in ABA bio- 
synthesis showed that VP14 (ZmNCEDl) and ZmA02 
(the homolog of Arabidopsis AA03) were significantiy up- 
regulated in CR and SR (Figure 4B). Several salt stress- 
inducible genes were identified in this study and showed 
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Figure 4 (See legend on next page.) 
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(See figure on previous page.) 

Figure 4 Differential gene expression of ZmPP genes and overall signaling pathway under salt stress. A Expression profiles of ZmPPs 
under salt stress in crown root (CR), primary root (PR), and seminal root (SR). Log2 signal intensities were used to create the heatmap. White 
mean only expressed in control or have no expression value both in control and salt stress treatment B Overview of salt-responsive systems in 
plants. C Real-time PGR results of salt-induced pathway components. Grey, control; black, salt treatment. 



in Additional file 13. To detect whether all of these genes 
were implicated in salt signal transduction pathways, 
RNA-seq analysis of three types of roots of maize under 
salt treatment was performed (Additional file 14). As the 
partner of PPs, MAPK was reported to function signifi- 
cantly in salt signaling in yeast [51]. To investigate the role 
of MAPK cascades in maize salt stress, the expression pat- 
terns of MAPK pathway components in the three types of 
roots were analyzed. As shown in Additional file 15, we 
found that ZmRaf46 was expressed at relatively high level 
in SR and CR compared with PR, indicating that MAPK 
cascades might play crucial roles in salt stress response in 
maize. Apart from MAPK, SNFl-related protein kinases 
(SnRKs) which acted as the downstream elements of PPs, 
function in salt stress tolerance as well [52]. Among the 14 
maize SnRK2s, ZmSnRK2.7 and ZmSnRK2.11 seem to play 
roles in salt stress tolerance with relative up-regulations 
(Additional file 14). Transcription factors (TPs) including 
MYB/MYC, NAC, WRKY, and zinc-finger protein have also 
been identified as salt stress-responsive factors in plants 
[50]. ZmWRKY33 (named ZmWRKY92.1 in our study), was 
reported to be significantiy up-regulated in salt stress [53,54]. 
In contrast, the expression level of ZmWRKY92.1 has 
no obvious change in CR under salt stress, but is 
slightly down-regulated in PR and SR, indicating that 
ZmWRKY92. 1 might act as a negative regulator of salt 
stress response. On the other hand, bZIP TPs function 
as ABA-dependent transcription factors and may be can- 
didate genes enhancing crops stress tolerance [33,55]. 
With the result of RNA-seq analysis, ZmbZIP23 expres- 
sion was up-regulated by salt treatment in maize CR and 
PR, suggesting that ZmbZIP23 might play a partial role in 
response to salt stress. MYB TPs ,which are characterized 
by highly conserved MYB DNA-binding domain, are in- 
volved in the regulation of salt stress response through 
ABA-independent pathway [50]. In this study, both 
ZmMYBOOS and ZmMYBlOO showed increased expressions 
in salinity, suggesting that they may act as master switch 
in the salt stress tolerance. NAC-type TPs also regulate 
some salt-responsive genes through ABA-dependent sig- 
naling pathway. In our analysis, high salinity stress induces 
several ZmNACs as well. It is notable that ZmNAC070.1 
and ZmNAC072 showed significantly up-regulated expres- 
sion in CR and SR. Additionally, plant exposure to high 
levels of NaCl also creates ionic stress in the form of cellu- 
lar accumulation of CI" and, in particular, Na^ ions. Previ- 
ous study showed that the vacuolar Na*/H* antiporter 
gene, ZmNHX displayed an up-regulation in root when 



exposed to high NaCl concentration [50]. What's more, 
salt stress induce either influx or efflux of nutrient ions 
such as Ca^^, K^, and NO^ . The salt overly sensitive 
(SOS) signaling pathway that comprises SOSl, SOS2, and 
SOS3 is a well known signaling pathway for resistance to 
salt stress [50]. SOS2, a Ser/Thr PK, forms a complex with 
SOS3, an EF-hand Ca ^-binding protein, to activate SOSl, 
which plays a central role in sodium extrusion and con- 
trolling long distance Na^ transport. Moreover, SOS2 also 
works to activate and interact with vacuolar N*/H^ and 
H*/Ca^^ antiporters and V-ATPase, which will contrib- 
ute to sequestration of excess Na* ions, further resulting 
in Na^ ion homeostasis. 

To confirm the expression patterns of these predicted 
genes under salt stress, 16 genes were selected for real- 
time PCR analysis. As shown in Figure 4C, most of the se- 
lected genes were up-regulated in salt stress, while few 
were down-regulated. Among these 16 genes, ZmMYBOOS, 
VP14, and ZmRaf46 were specially expressed under salt 
stress in both CR and SR, and ZmNAC070.1 displayed up- 
regulation in CR and had specific expression in SR. More- 
over, ZmS0S2 was expressed specially in CR but has no 
obvious fluctuation in SR. Taken together, real-time 
PCR-based expression profiling for these selected genes 
confirmed the outcome of RNA-seq analysis. Combined 
with RNA-seq, real-time PCR, and co-expression analysis 
(data not shown), ZmRaf46/49/17-ZmMKK3-ZmMPK3 
might represent a potential MAPK cascade involved in 
salinity stress. 

Differential expression of ZmPP genes under drought stress 

It is well known that drought is a major environmental 
factor determining plant productivity and distribution. 
The previous researches suggested that some PP genes 
not only played crucial roles in stress-induced ABA ac- 
cumulation but also function in redox signaling in maize 
[19,21]. PTP is a group of PP with 29 members in maize. 
To gain an insight into the functions of PTPs in ABA ac- 
cumulation and stress responses, we examined the effect 
of PAO (a PTPase-specific inhibitor) [21] on the expres- 
sion of the key enzymes involved in ABA biosynthesis, 
such as AA03, NCED3, SDRl, ABA3, and ZEP in Arabi- 
dopsis. ZmPPl is sensitive to the PAO and 50 ^iM PAO 
was use to inactivated the enzyme. We found that PAO 
can arrest the expression of genes encoding the key en- 
zymes {AA03, NCED3, and ABA3) in the ABA biosyn- 
thesis pathway, indicating that PTPs may have the 
potential roles in ABA biosynthesis (Figure 5A). Li et al. 
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(See figure on previous page.) 

Figure 5 Differential gene expression and real-time PCR analysis of ZmPP genes. A Effect of PAO on the key enzymes of ABA biosynthesis 
under water stress condition in Arabidopsis. 1, control; 2, PAO-treated; 3, negative control. B 1, stable expression of ZmPPl-GUS fusion protein in 
Arabidopsis chloroplast; 2, transient expression of ZmPPl-GFP fusion protein in onion epidermal cell; 3, transient expression of GFP protein only. 
C Analysis of gene expression in Arabidopsis transformed with estradiol-inducible ZmPPl. 1, wild-ethanol treated; 2, wild-estradiol treated; 3, 
ZmPPl over-expression 1 (ZmPP/Of /)-estradiol treated; ZmPP/Of2-estradiol treated. D Expression profiles of ZmPP genes under moderate drought 
stress (M/C), severe drought stress (S/C), and re-watering (R/C) as compared to control seedlings in Han21 and Ye478, respectively. E Real-time 
PCR analysis of 24 differently expressed PP genes under drought stress. F ABA biosynthesis and ABA-dependent pathway of response to drought 
in plants. The red arrows represent expression levels only detected by RNA-seq under drought stress. 



has showed that ZmRIPl (named ZmPPl in this study), 
which is a member of PTPs, functioned in redox signal- 
ing in maize [21]. To better understand the cellular 
function of ZmPPl, we further examined its subcellular 
location (Figure 5B, Additional file 16). When GFP- 
ZmPPl was introduced into onion epidermal cells, we 
found that ZmPPl had multiple subcellular locations, 
whereas GFP alone was localized to the nucleus and 
cytoplasm (Figure 5B). Moreover, histochemical staining 
revealed that GUS expression occurred mainly in the 
chloroplasts. However, stable expression of ZmPPl in 
Arabidopsis based on estrogen-receptor-based inducible 
system showed that ZmPPl is not a mediator of stress- 
induced ABA accumulation (Figure 5C). 

To test whether the ZmPTPs play roles in drought 
stress, microarray data of different drought treatments 
were used to generate a heatmap. This genomic study of 
maize gene expression in response to water deficit include 
two different inbred lines, namely, drought-tolerant line 
Han21 and drought-sensitive line Ye478 (Figure 5D). Be- 
sides DNA microarray analysis, RNA-seq analysis was per- 
formed to investigate the functions of ZmPTPs in different 
drought sensitive organs, namely, leaf and cob, under 
drought and well-watered conditions [56]. Imposing the 
two-fold change requirement with a P value <0.05, the 
public microarray data showed that expression of ZmPP67, 
a member of PTP subfamily, was up-regulated in Han21. 
Further observation from RNA-seq data showed that 
ZmPP67 was significantly up-regulated in both leaf and 
cob, indicative of its important role in response to water 
deficit conditions. To confirm the expression patterns of 
ZmPTPs under water deficit condition, real-time PCR ana- 
lysis was carried out using 8 differently expressed genes. 
As shown in Figure 5E, ZmPP67 was obviously up- 
regulated in water deficient, which is consistent with our 
previous assumption that ZmPP67 might have an essential 
part during drought stress. Moreover, from RNA-seq and 
real-time PCR data, it was interesting to find that the other 
two PTP genes, ZmPP54 and ZmPP59, whose expression 
levels were also up-regulated. Of these 8 differentially 
expressed genes, 4 genes {ZmPP12S, -133, -101 and -68) 
were down-regulated in water deficit condition. But the 
functions of these 4 genes in drought stress are still 
unlcnown. 



Previous studies have reported that 6 Arabidopsis Group 
A PP2Cs interacted with ABA receptors and functioned as 
key negative regulators in ABA signaling [9,57]. Thus, in 
addition to PTPs, other PP genes may be involved in drought 
response. In RNA-seq data (Additional file 17), 10 genes 
{ZmPP6, -21, -29, -31, -113, -124, -127, -130, -146, 
and -154) were significantly up-regulated in leaf with 
the criterion that fold change > 2. Of these 10 genes, all 
of them (except ZmPP124) were belong to PP2C class. 
After integrating the RNA-seq and real-time PCR data 
(Figure 5E, Additional file 17), it is notable that several 
putative PP orthologs among Arabidopsis, rice and 
maize showed strikingly consistent expression patterns, 
which lend further supporting to the existence of func- 
tional conservation among these species. For instance, 
ZmPP2C (named ZmPP76), which was reported to be 
essential for plant stress signal transduction, showing 
up-regulated expression pattern in drought stress [19]. 
In our study, ZmPP76 was slightly up-regulated, whose 
expression pattern was in good agreement with the pre- 
vious report [19]. Interestingly, another stress-induced 
gene ZmPP2Ca, named ZmPP8 in this study, was up- 
regulated by drought. In contrast, previous study has re- 
vealed that down-regulation of the ZmPP2Ca transcript 
was observed in drought tolerant lines, whereas tran- 
script levels were up-regulated in sensitive lines under 
drought stress [17]. In addition, ZmPPl 12, the homolog 
of rice protein phosphatase OsPP2Cl, was up-regulated 
more significantly in the drought sensitive line Ye478 
than that in the drought tolerant line Han21. Previous 
study showed that OsPP2Cl was highly induced in 
response to abiotic stresses, even up-regulated at 
low temperature and drought conditions [58]. On the 
contrary, RNA-seq data and real-time PCR analyses of 
ZmPPl 12 exhibited decreased expression level in drought. 
Recently, Americo et al. reported that two members of 
Group A PP2C (ABIl and PP2CA) can interact with the 
SnRKl catalytic subunit, causing its dephosphorylation 
and inactivation, while the inhibition of PP2C allows ABA 
to promote SnRKl activation [59]. ABI2, the homolog of 
ABIl in Arabidopsis, has been reported to involve in ABA 
signal transduction with up-regulated expression pattern 
[60]. Similarly, ZmPP49, the homolog of Arabidopsis ABI2 
in maize, was up-regulated, indicating that ZmPP49 might 
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interact with ABA receptors to suppress the activity of 
SnRKs. 

Apart from the above putative genes up-regulated by 
ABA, another two genes {ZmPP138 and ZmPP82) were 
expressed at nearly equivalent levels. As previously re- 
ported, the expression level of ZmPP138 (published 
ZMPP2) was relatively high in root and etiolated shoots 
[20]. Also, OsBIPP2Cl, the homolog of ZmPP82 in rice, 
was up-regulated upon drought treatments [61]. In an 
attempt to examine the effects of drought stress on the 
expression of the PP2A genes, comprehensive analysis 
on microarray, RNA-seq and real-time PCR was done 
(Figure 5D-E, Additional file 17). Especially, ZmPP17 
and ZmPP79, two homologs of rice PP2As {OsPP2A-4 
and OsPP2A-2), showed special expression patterns. 
OsPP2A-4 is expressed at comparably high levels in 
stems, roots and flowers, while an appreciably lower level 
of expression is observed in leaves. However, ZmPPl 7 was 
not expressed during drought in leaf [62]. Additionally, 
the expression level of OsPP2A-2 was significantiy down- 
regulated in leaves and stems in response to drought 
stress, which is consistent with the expression pattern of 
ZmPP79 in our study [62]. Another PP2A gene, ZmPP6S, 
showed down-regulated expression in leaf under drought 
stress. In addition to those genes that have been reported 
or had homologs in other species, some maize PP genes 
were also found to be up-regulated under water deficit, 
such as ZmPP29, -113, -146 and -154. Typically, 
ZmPP154 is a member of Group A PP2C in maize. In 
addition, ZmPP7, ZmPPlOS and ZmPP40 were especially 
down-regulated under water deficit condition. 

A putative schematic of the ABA biosynthesis and 
ABA cellular response during drought stress was pro- 
posed based on the described web-like networks of ABA 
biosynthesis and ABA signaling under water shortage 
condition (Figure 5F). Compared with Arabidopsis, 
maize has less characterized drought-induced genes 
(Additional files 18 and 19). Based on our investigation 
data and the previous research results [63,64], we found 
that ABA biosynthetase genes were up-regulated in the 
drought- tolerant line Han21, although their protein 
levels have not been examined in most cases. 

Expression profiling of ZmPP genes under cold stress 

Cold responses have been observed in many plants, like 
Arabidopsis and rice, which alter mRNA levels of genes 
belonging to multiple independent pathways to response 
to low temperature [65,66]. To search for the molecular 
mechanism of high intraspecific diversity in cold tolerance 
of maize, transcriptomic analysis was studied in two inbred 
lines contrasting in chilling tolerance (Additional file 20). 
With the criterion that the response must be at least 
2.8-fold (|log2[(expression level in the cold):(expression 
level in control conditions)]] > 1.5) [11], 13 ZmPP genes 



{ZmPP6, -24, -29, -66, -77, -82, -92, -112, -116, 
-127, -149, -154, and -155) were found to be significandy 
up-regulated, while three genes {ZmPP69, -87, and -101) 
displayed obvious down-regulation. AtPP2CA, an Arabi- 
dopsis group A PP2C, was reported as a negative regulator 
in ABA responses during cold acclimation [12]. ZmPP6, 
the homolog of AtPP2CA in maize, which was up- 
regulated in cold-tolerant line ETH-DH7. In addition, 
ZmPP2C2, named ZmPP159 in this study, occurred to be a 
positive regulator of cold resistance in plant [18], whose ex- 
pression was slightly up-regulated. It is conceivable that 
maize PP2C activity could be modulated via cold stress- 
induced response. To investigate the role of MAPK cas- 
cades in cold tolerance, we performed microarray analysis 
of components of MAPK pathway (Additional file 21). In 
cold tolerant line, three MAPKs (ZmMPK6, -8 and -12), 
two MAPKKs {ZmMKK2 and ZmMKK4) and four MP3Ks 
{ZmMAPKKKl, ZmMAPKKK17, ZmZIK9 and ZmRaf43) 
were significantly up-regulated. To explore the gene ex- 
pression patterns in our maize inbred lines under cold 
stress, the obviously up-regulated ZmPPs and the compo- 
nents of MAPK cascades were selected for real-time PCR 
analysis. As shown in Additional file 22, all the selected 
genes were up-regulated under cold treatment. We pro- 
posed that MAP3Kl/17/ZIK9/Raf43-MKK2/4-MAPK6/8/ 
12 might represent putative MAPKK-MAPK or MAPKKK- 
MAPKK-MAPK cascades based on microarray, real-time 
PCR and co-expression analyses (data not shown). The 
predicted MAPKKK-MAPKK-MAPK phosphorylation 
regulatory network constitutes a valuable resource to 
understand the function and specificity of MAPK signaling 
systems in maize. It is worth mentioning that ZniMP3Kl- 
and ZmMPK12-overexpression in tobacco showed higher 
seed germination rate under low temperature treatment 
(12°C) (data not shown). 

Promoter analysis and miRNA targets ZmPPs 

In genetics, a promoter is a region of DNA that initiates 
transcription of a particular gene and confers developmen- 
tal and/or environmental regulation of gene expression. 
Discovery of cw-regulatory elements in the promoter re- 
gions is essential to understand the spatial and temporal 
expression pattern of ZmPP genes. To further expound 
the transcription regulation machineries of the 159 ZmPP 
genes in maize, a comprehensive promoter analysis was 
performed by PlantCARE database [67]. The occurrence 
of c«-elements in ZmPP genes are showed in Additional 
file 23. The statistical analysis showed that 12 cw-elements 
were enriched in maize PP genes, which were involved in 
regulation of gene expression under stress conditions. 
Nearly all genes (except ZmPP13, -71, -125) contained 
part of these 12 ci's-elements, and the elements were vari- 
ably present in the promoter regions of most ZmPP genes. 
Interestingly, we found that the majority of the ZmPP 
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genes contained ABRE, G-Box and MBS, which are the TF 
binding sites for bZIP, GBF (G-Box-binding factor) and 
MYB, respectively. The data might indicate ds-elements 
in this study may play major roles in regulating the expres- 
sion of PP genes in response to different stresses in maize. 

MicroRNAs (miRNAs) play important roles in plant 
post-transcriptional gene regulation by either cleaving 
mRNA transcripts or repressing protein translation [68]. 
Putative miRNAs targeting the ZmPP genes were identi- 
fied using Target-align program. It showed that about 35 
ZmPP genes were targeted by maize miRNAs (Additional 
file 24). In addition to regulating growth and development, 
most miRNAs in plants play important roles in the regula- 
tion of various cellular processes underlying plant adapta- 
tion to environmental stresses. Among these target genes, 
6 and 8 members belong to Group A PP2C and PP2A, re- 
spectively, indicating that miRNAs in the present study 
would assist in understanding the post-transcriptional 
control of gene regulation during physiological and stress- 
induced cellular responses. 



Conclusions 

Conclusively, our results present a comprehensive account 
of PP-encoding genes and provide new insights into the 
phylogenetic relationships and characteristic functions of 
maize PPs. In addition, analyses of expression profiles 
based on microarray and RNA-seq method unraveled 
their probable functions during different developmental 
stages and abiotic stresses and will be useful in studies 
aimed at revealing the global regulatory network in maize 
development and abiotic stress responses, thereby contrib- 
uting to the maize molecular breeding with enhanced 
quality traits. 

Additional files 



Additional file 1: Figure 51. Phylogenetic analysis of maize, rice, and 
Arabidopsis PP genes. An un-rooted NJ tree is made based on the 
catalytic domain sequences of maize, rice, and Arabidopsis PPs. PPs from 
Arabidopsis, rice, and maize belong to the same class falling in the same 
clades. Scale bar represents 0.1 amino acid substitutions per site. 

Additional file 2: Figure 52. The map of intron/exon arrangement of 
ZmPP genes. 

Additional file 3: Figure 53. Qrcos diagram of protein phosphatase 
gene pairs between maize and rice genomes. Outer two circles 
Distribution of each of the PP genes and scaled chromosomes for each 
species in million bp (Mb) units, respectively. Histograms below each 
chromosome Number of introns of PK genes: green < 5 introns, red > 10 
introns. Boxes Syntenic regions. Colors are assigned to the syntenic 
regions according to the colors of the corresponding chromosome. 
Innermost colored lines interconnect putative orthologous PP gene pairs 
between rice and maize. 

Additional file 4: Figure 54. Heatmap showing the clustering of ZmPPs 
according to their expression profiles of 60 detected transcripts at 
different stages/organs of maize. Red, white and green indicate high, 
medium and low levels of gene expression, respectively. E2enzyme was 
used as a internal control. 



Additional file 5: Figure 55. A putative schematic of root development 
signaling pathway in maize. The little colored blocks besides the gene 
expression level of pathway components under different developmental 
stages. 

Additional file 6: Figure 56. RNA-seq differentially expressed analysis 
of ZmPP genes along maize leaf developmental gradients, namely, 
base, -1 cm, +4 cm, and tip. E2enzyme was used as a internal control. 

Additional file 7: Figure 57. Expression patterns of ZmPP genes in 12 
diverse maize reproductive tissues. DAP, days after pollination; Pre-em, 
preemergence; Post-em, postemergence. E2enzyme was used as a 
internal control. 

Additional file 8: Figure 58. Expression profiles of signaling 
components under salt stress in CR, PR, and SR, respectively. Log2 based 
fold changes were used to create the heatmap. Expression values 
highlight with white mean only expressed in control or have no 
expression value both in control and salt stress treatment 

Additional file 9: Figure 59. The vector map of pEZS-NL-ZmPP/. 

Additional file 10: Figure 510. Real-time PCR analysis of representative 
ZmPPs and MAPK-cascade genes under cold treatments. The maize 
ACTIN? gene was used as endogenous control to normalize data. 
E2 enzyme was used as a internal control. 

Additional file 11: Table 51. List of primers used in this study. 

Additional file 12: Table 52. The identified maize protein 
phosphatases and their related information. 

Additional file 13: Table 53. List of 152 detected ZmPP genes and 
their expression patterns at 18 selected tissues. 

Additional file 14: Table 54. List of putative root development 
signaling components in maize. 

Additional file 1 5: Table 55. List of flowering time signaling 
components in maize. 

Additional file 16: Table 56. List of putative salt stress signaling 
components in maize. 

Additional file 17: Table 57. List of expression values of MARK 
cascade-related genes and key genes of ABA biosynthesis in maize under 
salt stress. 

Additional file 18: Table 58. List of FPKM values of ZmPPs in maize 
reproductive (cob) and vegetative tissue (leaf) under both drought and 
well-watered conditions. MCC and MCD stand for maize ovary, well 
watered and drought, respectively; MLC and MLD stand for maize basal 
leaf meristem, well watered and drought, respectively. Numbers 1 and 2 
indicate the two biological replicates. The extent of differential expression 
is measured in terms of fold change and (-) indicates failure to calculate 
or undetected values. Values in red and blue indicate the fold increase 
and decrease in expression in the drought-stressed tissue, respectively. 

Additional file 19: Table 59. List of EPKM values of components involved 
in ABA biosynthesis and ABA-dependent pathway of response to drought in 
vegetative tissue (leaf) under both drought and well-watered conditions. 
MLC and MLD stand for maize basal leaf meristem, well watered and 
drought respectively. Numbers 1 and 2 indicate the two biological 
replicates. The extent of differential expression is measured in terms of fold 
change and (-) indicates failure to calculate or undetected values. Values in 
red and blue indicate the fold increase and decrease in expression in the 
drought-stressed tissue, respectively. 

Additional file 20: Table 510. List of putative components involved in 
ABA biosynthesis and ABA-dependent pathway of response to drought 
in maize. 

Additional file 21: Table 511. List of expression values of ZmPP genes 
under cold stress. ETH-DH7 log(c), ETH-DH7 log(k), ETH-DL3 log(c), 
ETH-DL3 log(k) means intensity (AU) of fluorescence of labeled aaRNA 
from cold-treatment (c) or control (k) hybridizing to the probe, shown as 
log2 in the respective maize line. DH7 log(c/k), DL3 log(c/k) means log2 
of ratio of expression in cold-treated vs. control of the probe in the 
respective maize line, respectively. Values in red and blue indicate the 
fold increase and decrease in expression in the drought-stressed tissue, 
respectively. 
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Additional file 22: Table SI 2. List of expression values of IVIAPK 
cascades genes in cold stress. ETH-DH7 log(c), ETH-DH7 log{k), ETH-DL3 
log(c), ETH-DL3 log(k) means intensity (AU) of fluorescence of labeled 
aaRNA from cold-treatment (c) or control (k) hybridizing to the probe, 
shown as log2 in the respective maize line. ETH-DH7 log(c/k), ETH-DL3 
log(c/k) means log2 of ratio of expression in cold-treated vs. control of 
the probe in the respective maize line, respectively. 

Additional file 23: Table SI 3. List of promoters in ZmPP genes. 

Additional file 24: Table S14. List of predicted mlRNA-regulated ZmPP 
genes. 
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